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Using an efficient numerical scheme that exploits spatial symmetries and spin- 
parity, we have obtained the exact low-lying eigenstates of exchange Hamihonians 
for ferric wheels up to Fei2. The largest calculation involves the Fei2 ring which 
spans a Hilbert space dimension of about 145 million for Ms=0 subspace. Our cal- 
culated gaps from the singlet ground state to the excited triplet state agrees well 
with the experimentally measured values. Study of the static structure factor shows 
that the ground state is spontaneously dimerized for ferric wheels. Spin states of 
ferric wheels can be viewed as quantized states of a rigid rotor with the gap between 
the ground and the first excited state defining the inverse of moment of inertia. We 
have studied the quantum dynamics of Feio as a representative of ferric wheels. We 
use the low-lying states of Feio to solve exactly the time-dependent Schrodinger 
equation and find the magnetization of the molecule in the presence of an alter- 
nating magnetic field at zero temperature. We observe a nontrivial oscillation of 
magnetization which is dependent on the amplitude of the ac field. We have also 
studied the torque response of Fei2 as a function of magnetic field, which clearly 
shows spin-state crossover. 

PACS numbers: 75.45. -t-j, 75.50.Xx, 75.60.Ej 



I. INTRODUCTION 

In recent years polyoxometalates which are practical realization of nanomagnets have become the 
area of intense research because of their enormous variety of structure and fascinating magnetic 
properties. A particular aesthetic class is that of the ring-shaped iron(III) compounds denoted 
as molecular ferric wheels. The decanuclear wheel Feio synthesized by Lippard et. al.0 may be 
regarded as a prototype of this class. Meanwhile synthesis of ferric wheels with 6, 8, 12 and 18 
have also been reportedl. These materials have dominant anti- ferromagnetic coupling between 
Fe(III) spins and a singlet ground state. The magnetic properties of such nanoscopic molecules 
result from the interplay of superexchange interactions between the atomic spins, dipolar coupling 
of the local moments and on-site anisotropics arising from ligand configurations. The emergence 
of ferric wheels has led to a renewed interest in the properties of the Heisenberg chain, especially 
for large spin values. In 1964 Fisher used a classical treatment^, in which spin quantization 
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is absent, to study the properties of Heisenberg Hamiltonian and found analytical solutions for 
the thermodynamic properties of the system. Numerous quantum mechanical calculations have 
been made for Heisenberg Hamiltonian, mostly for spin- 1/2 chains. Calculations for larger S 
became more interesting after Haldane's conjectureS regarding the presence of an energy gap in 
the excitation spectrum from the ground state for integer S chains. 

To fit the experimental temperature dependence of the magnetic susceptibility data for Feio Lip- 
pard et. al. adopted a classical spin treatment to obtain the value of exchange interaction parameter 
J. But below 50 K this treatment fails. Ferric wheels with S=5/2, and system size of up to 8 sites 
have been treated exactly using irreducible tensor operator techniquJ^ with the aid of point group 
symmetry as well as by using the invariance of spin Hamiltonian with regard to interchange of 
spin sitesi. Recently there have been theoretical studies to explain the low temperature magnetic 
susceptibility data of FciqE The magnetization of ferric wheels exhibit a step like field dependence 
at low temperature due to the occurrence of field induced ground state level crossing, a spectacular 
manifestation of quantum size effect in these nanomagnetsE'i. There have been NMR and specific 
heat studies of these ferric wheels to investigate the energy level structurc§0. In appropriate 
parameter regime these ring systems are considered to be the candidates for the observation of 
macroscopic quantum phenomena, in the form of quantum coherent tunneling of the Neel vectoiO. 
To understand these low temperature properties of ferric wheels in detail we need to know the low 
lying eigenvalue spectrum for these systems. It is also of interest to compare and contrast the zero 
temperature density of state of S=5/2 rings with that of the exactly known S=l/2 chain. In this 
paper we have used spin parity together with rotational symmetry of the ferric wheels to obtain 
the low-lying eigenvalue spectrum for rings with 6, 8, 10 and 12 spin-5/2. The dispersion spectrum 
reveals interesting features. There have been recent reports of solving the low-lying eigen-spectrum 
of Feio using the Density Matrix Renormalization GroupEl. To the best of our knowledge ours is 
the first study of of ferric wheels using exact methods up to a ring size of 12, S=5/2 iron(III) ions. 

We have also studied the quantum spin dynamics of ferric wheels in presence of an alternating 
field after setting up the Hamiltonian matrix in the subspace of low-lying states. This Hamiltonian 
includes multipolar spin-spin interaction terms besides a time varying magnetic field (Zeeman 
term). We have then evolved an initial state, which is taken to be the ground state with a specific 
value of Ms (the z-component of the total spin) in the absence of the magnetic field, by using the 
time-dependent formulation of the problem in the restricted subspace. We observe a nontrivial 
oscillation of the magnetization whose frequency depends on the amplitude of the alternating field. 
This phenomena is very similar to what has been observed in case of uniaxial magnets0. 
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II. MODEL HAMILTONIAN AND COMPUTATIONAL DETAILS 



A. Symmetry Adaptation of Correlated States 



The full symmetry of electronic states is of central interest in quantum theory. Since nonrela- 
tivistic Hamiltonians H does not depend explicitly on spin, molecular eigenstates can be labeled 
by the total spin S and the appropriate irreducible representation of the point group. The model 
Hamiltonian employed in this study is the isotropic exchange Hamiltonian involving exchange 
interaction between nearest neighbors, 



where the exchange interaction takes the values dictated by experimental studies of structure 
and magnetic properties. For Feg the J value depends on the central alkali metal atom : for NaiFce 
J=32.77 K whereas for LiiFeg J=20.83 K. In Fes, Feio and Fei2 J is 22 K,15.56 K and 31.97 K 
respectivelyEi. 

The total dimensionality of the Fock space of the ferric wheel is given by 



where 'n' is the total number of spins in the wheel and Si is the spin on each ion. In case of Fcio, 
there are 10 iron(III) ions with the dimensionality of the Fock space 60,466,176. In case of Fei2, 
which we have solved exactly this dimension is 2,176,782,336. 

Specializing to given total Ms leads to Hilbert space dimensionalities which are lower than the 
Fock space dimensionality. In the case of the Fei2 cluster the Ms = space has a dimensionality 
of about 145 million (144,840,476). The major challenge in exact computation of the eigenvalues, 
and properties of these spin clusters lies in handling such large bases and the associated matrices. 
While the dimensions look overwhelming, the matrices that represent the operators in these spaces 
are rather sparse. Usually, the number of nonzero elements in a row is of the order of the number 
of exchange constants in the Hamiltonian. This sparseness of the matrices allows one to handle 
fairly large systems. However, in the case of spin problems, generating the bases states and using 
the symmetries of the problem is nontrivial. 

The isotropic exchange Hamiltonians conserve total spin, S, besides the z-component of the total 
spin, Ms. Besides these symmetries, the geometry of the cluster also leads to spatial symmetries 
which can often be exploited. The simplest way of generating bases functions which conserve total 
spin is the VB method that employs the Rumer-Pauling rulJll. It is quite easy to generalize the 
Rumer-Pauling rules to a cluster consisting of objects with different spins to obtain states with 
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desired total spin, S. However, setting-up the Hamiltonian matrix in such a basis can be computa- 
tionaUy intensive since the exchange operators operating on a "legal" VB diagram (diagram that 
obeys Rumer-Pauling rules) could lead to "illegal" VB diagrams and resolving these "illegal" VB 
diagrams into "legal" diagrams would present the major bottle-neck. Indeed, the same difficulty is 
encountered when spatial symmetry operators operate on a VB function. Thus, the extended VB 
methods are not favored whenever one wishes to apply it to a motley collection of spins or when 
one wishes to exploit some general spatial symmetries that may exist in the cluster0. 

It is advantageous to partition the spaces into different total spin spaces because of the usually 
small energy gaps between total spin states which differ in S by unity. To avoid the difficulties 
involved in working with total spin eigenfunctions, we exploit parity symmetry in the systems. 
The parity operation involves changing the z-component of all the spins in the cluster from Mg; 

s s- 

to — Msj. There is an associated phase factor with this operation given by (-1)^'°'+Z^.=". The 
isotropic exchange operator remains invariant under this operation. If this symmetry is employed in 
the Ms = subspace, the subspace is divided into "even" and "odd" parity spaces depending upon 
the sign of the character under the irreducible representation of the parity group. The space which 
corresponds to even (odd) total spin we call the even (odd) parity space. Thus, employing parity 
allows partial spin symmetry adaptation which separates successive total spin spaces, without 
introducing the complications encountered in the VB bases. However, the VB method can lead to 
complete factorization of the spin space leading to smaller complete subspaces. 

In the ferric wheels, besides spin symmetries, there also exists spatial symmetries. The topology 
of the exchange interaction leads to a Cn point group symmetry, where 'n' is the number of iron 
ions in the ring. Hence, Fee, Feg, Feio, Fei2 will have Ce, Cg, Cio and C12 symmetry respectively. It 
should be mentioned that among these ferric wheels only Feio is strictly having a ten fold rotational 
symmetry, rest of them approximately have the above mentioned symmetry. For computational 
advantage we have assumed the rotational symmetry for all the ferric wheels. This point group 
appears at first site to present difficulties because the characters in the irreducible representation 
are in some cases complex which could lead to complex bases functions. This, however, can be 
avoided by recognizing that in the Cn group, states with wavevectors 'k' and '-k' are degenerate. 
We can therefore construct a linear combination of the 'k' and '-k' states which is real. The 
symmetry representations in the Cn group would then correspond to the labels A, B and E, 
with the characters in the E representation given by 2cos(rfc) under the symmetry operation CJ^, 
with k = 7r/n. The parity operation commutes with the spatial symmetry operations and the full 
point group of the system would then correspond to the direct product of the two groups. Since 
both parity and spatial symmetries can be easily incorporated in a constant Ms basis, we do not 
encounter the difhculties endemic to the VB theory. 
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The generation of the eomplete basis in a given Hilbert space requires a simple representation 
of a state on the computer. This is achieved by associating with every state a unique integer. 
In this integer, we associate n; bits with spin Si, such that nj is the smallest integer for which 
2"' > (2si + 1). In the integer that represents the state of the cluster, we ensure that these ni 
bits do not take values which lead to the ni bit integer value exceeding (2si + 1). For each of 
the allowed bit states of the Uj bit integer, we associate an value between — Si and Si. For a 

n 

spin cluster of 'n' spins, we scan all integers of bit length N = S ni and verify if it represents a 

i— 1 

basis state with the desired Ms value. In Fig. 1, we show a few basis functions with spcKiified 
Ms value for some typical ferric wheels along with their bit representations and the corresponding 
integers. Generation of the bases states is usually a very fast step, computationally. Generating 
the basis as an ordered sequence of integers that represent them also allows for rapid generation 
of the Hamiltonian matrix elements as will be seen later. 

Symmetrization of the basis by incorporating parity and spatial symmetries involves operation on 
the constant Mg = basis by the symmetry operators. Since spatial symmetry operators permute 
the positions of equivalent spins, every spatial symmetry operator operating on a basis function 
generates another basis function. Every symmetry operator can be represented by a correspondence 
vector whose i*'' entry gives the state that results from operating on the i*^ state by the chosen 
operator. This is also true for the parity operator, in the Ms = subspace. 

The first step in constructing symmetry adapted linear combinations is to represent the symmetry 
operators in the chosen basis as matrices. In our case, the symmetry operators are such that 
symmetry operation by any operator on a basis state leads to a resultant which is a single basis 
state. Thus all our symmetry operators can be represented as vectors; the entry in position 'i' 
gives the index of the basis function generated by symmetry operation on the basis state 'i'. Since 
the basis is very large, it is prohibitive to store and manipulate the full basis together with all 
the associated symmetry vectors. To avoid these difficulties, we have constructed the symmetry 
matrices in small invariant subspaces. These invariant subspaces are obtained by sequentially 
choosing a state and operating on it by all the symmetry operators. This gives rise to a set of 
states on which we again operate by all the symmetry operators, and continue this process until 
no new basis states are generated. The collection of all these basis states resulting at the end of 
this process is the invariant subspace. We can set-up symmetry combinations of the basis states in 
each of the invariant subspaces independently. The symmetry combinations can now be obtained 
operating on each state by the group theoretic projection operator, 

^r. = l^XvmR (3) 
on each of the basis states of the invariant subspace. Here Fj is the i*'' irreducible representation. 
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R is the symmetry operation of the group and xr; {R) is the character under R in the irreducible 
representation Ti . This process is repeated with the next basis state that has not appeared in any 
of the invariant subspaces aheady constructed. The process comes to an end when aU the basis 
states have appeared in any one of the invariant subspaces. 

The resuhing symmetrized basis is usuahy over complete in each of the invariant subspaces. 
The linear dependencies can be eliminated by a Gram-Schmidt orthonormalization procedure. 
However, in most cases, ensuring that a given basis function does not appear more than once 
in a symmetrized basis is sufficient to guarantee linear independence and weed out the linearly 
dependent states. A good check on the procedure is to ensure that the dimensionality of the 
symmetrized space in the invariant subspace agrees with that calculated from the traces of the 
reducible representation obtained from the matrices corresponding to the symmetry operators in 
the chosen invariant subspace. Besides, the sum of the dimensionalities of the symmetrized spaces 
should correspond to the dimensionality of the unsymmetrized invariant subspace in each of these 
subspaces. 

Generation of the Hamiltonian matrix is rather straight forward and involves operation of the 
Hamiltonian operator on the symmetry adapted basis. This results in the matrix SH, where S is 
the symmetrization matrix representing the operator Pr- and H is the matrix whose elements hij 
are defined by 

H|i>=Ehij|j > . (4) 
j 

The states {i} correspond to the unsymmetrized basis functions. The Hamiltonian matrix in the 
symmetrized basis is obtained by right multiplying the matrix SH by . The resulting symmetric 
Hamiltonian matrix is stored in the sparse matrix form and the matrix eigenvalue problem is solved 
using the Davidson algorithm0. 

Computation of the properties is easily done by transforming the eigenstate in the symmetrized 
basis into that in the unsymmetrized basis. Since the operation by any combination of spin 
operators on the unsymmetrized basis can be carried out, all relevant static properties in different 
eigenstates can be obtained in a straightforward manner. 

B. Quantum Spin Dynamics 

We have studied the dynamics in ferric wheels by setting up the Hamiltonian matrix in the 
desired Mg space, which in all cases is restricted to Ms = 0, 1 and 2. In each subspace we have 
obtained a few low-lying states using the Davidson algor ithmQ. We have also calculated the spin- 
spin correlation functions in each of the states. Using the spin-spin correlation functions, we have 
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computed the expectation value of S^^^^j operator, from which we have identified the total spin 
of the state. We observe from the eigenvalue spectrum of all the ferric wheels that ground state 
and the first excited state are spin singlet (S=0) and triplet (S=l) state respectively, they belong 
to two different spatial symmetry subspaces. So they will not mix unless there is a perturbation 
which spoils the Cio symmetry of the molecule. We also notice that the first (triplet) and the 
second (quintet) excited states again fall into different symmetry subspaces. This is true in all the 
ferric wheels we have studied. 

To study quantum dynamics we have considered the following HamiltonianS, 

3z, total (5) 



H — Eg — D S^^total + C (Sx, total ~ total) ~ g h(t) Sz 



Here D is the quadratic anisotropy factor, g is the Lande g-factors for the iron(III) spin respectively, 
and h(t) is the time-dependent magnetic field, expressed as h(t) — Hocos(wt), where uj is the 
frequency at which field is ramped and Hq is the amplitude of the field. We have chosen D = 
8.8 X 10^'^ and c = 10^"^ (in units of J) in accordance with the experimental values for Fciq. We 
take g — 2.0. The constants Eg in Eq. (||) correspond to the lowest energies obtained from Eq. 
(H). The second order anisotropy term allows transitions between states with AMs = ±2. Both 
second and third terms in Eq. (|^) arises due to magnetocrystalline anisotropy. The exact form of 
the anisotropy in ferric wheel is not very well established. We have included anisotropy terms only 
up to second order in the spin variables. The anisotropy in the plane can be formed artificially, 
e.g. by means of external electric or magnetic fields, pressure or using anisotropic substrateEl. 

The Hamiltonian in Eq. (|^) does not allow the mixing of the (singlet) ground state, (triplet) 
first excited state and (quintet) second excited state in Feio because of symmetry reason already 
mentioned earlier. To observe spin dynamics we apply a constant magnetic field that leads to 
triplet being the ground state of the system. This constant shifting field plays here the role of 
a chemical potential and is used to provide the optimal conditions for observing the considered 
quantum spin dynamics. We can represent the external field as 

B = Bo + Hocos(tjt) (6) 

To study the evolution of magnetization as a function of applied oscillatory field we start from a 
fixed axial field amplitude (in units of 3/h) and time is increased in small time step of At =0.01 
(/i/J), frequency of the field is kept constant. 

We have studied the time evolution of the system by solving the time-dependent Schrodinger 
equation, 

i;.^ = H(t)v . (7) 
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We assume the system to be in the state with = 1, Ms = — 1 at time t = 0. This is the initial 
state which is time evolved according to the equation 

^{t + At) = e-»^(*+¥)^*A ^(i) . (8) 

The evolution is carried out by explicit diagonalization of the Hamiltonian matrix H{t + ■^), 
and using the resulting eigenvalues and eigenvectors to evaluate the matrix of the time evolution 
operator e^*^(*+^''^*/''. We set up the Hamiltonian matrix for time evolution in the truncated 
basis of 3 states corresponding to total spin 5=1. We repeatedly carry out the the time evolution 
in small time steps of size At to obtain time evolution over longer periods. 



III. RESULTS AND DISCUSSION 
A. Analysis of Low-lying Spectrum 

We have solved the exchange Hamiltonian (Eq. (|l])) exactly using the method mentioned earlier 
to obtain the low-lying eigenvalue spectrum for 6, 8, 10 and 12 site iron(III) rings. We find that 
ground state, first, second and third excited states are respectively singlet, triplet, quintet and 
heptet for all ferric wheels. We notice that there is no accidental degeneracy between the energy 
levels belonging to different symmetry subspaces. Ground state switches between A+ (k =0) and 
(k = tt) subspace for even and odd value of 'm' respectively, in N=2m ring systems. This was 
also observed by Mattheiss in case of a spin-1/2 chain and can be understood from Marshall's sign 
rule. The gap between the ground state and the first excited state is shown in Fig. 2 as a function 
of inverse ring size. According to Haldane conjecture, the gap should extrapolate to zero. The 
extrapolated value, while small is still finite, suggesting that in these rings the finite size effects 
are still at play at the ring sizes we have studied. 

Using the exchange constants estimated for the different ring systems, we estimate the gap 
between the ground state and the first excited state to be 22.67 K, 11.81 K, and 6.88 K for Na : Feg, 
Fes, Feio respectively. Our calculated values compare very well with the experimental values, which 
are 22.0 K, 12.1 K and 6.45 K for Na : Feg, Feg and Feio respectivelyQ. This agreement shows that 
for all practical purposes ferric wheels can safely be assumed as rings neglecting the slight deviation 
from the exact circular geometry. Our calculated gap for and Fei2 is 12.09 K corresponding to 
the exchange constants predicted from experiments. However, experimental estimate of this gap 
is lacking. 

If we define Si to be the energy difference between the i-th excited state and the ground state. 
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we find from Fig. 3 that the following relationship is satisfied for the ferric wheels, 

= X E, (9) 

where Ei is the energy gap between the ground state and the first excited state. This indicates 
that the lowest spin state obeys the Lande interval rule, in agreement with Taft's conjecture. If 
we assume Ei to be the inverse of the moment of inertia, then the above expression gives the 
rotational energy of a rigid rotor in a state with quantum number Si. Thus, the spin states of 
ferric wheels can be viewed as the quantized states of a rigid rotor. 

In Fig. 4 we have shown the dispersion spectrum of ferric wheels. The value of k corresponding 
to a wavefunction can be defined as, 

T^p = e'V (10) 

where T is the translation operator which in case of a ring rotates the ring by one lattice spacing. 
We have used the spatial C„ symmetry of the ferric wheels, which enables us to identify the value 
of k easily for a specified eigenfunction. We find that, in the momentum (fc) sector which contains 
the ground state, the lowest excitation is to a quintet state. For all the other fc-sectors the lowest 
excitation is to a triplet state. Previous studieJll on antiferromagnetic spin-1/2 Heisenberg chains 
show that the excitation spectrum is given by fku — (7r/2) J| sinfcj, where k is the wave vector of 
the excited states measured with respect to that of the ground state. The simple antiferromagnetic 
spin-wave theory based on the use of the Holstein-Primakoff transformation for each sublattice, 
leads to excitation spectrum {S is the magnitude of the individual spin Si), 

hu = 2 J S \ smk\. (11) 

This relation is supposed to be correct for — > oo. We notice that the excitation spectrum for 
ferric wheels can be fitted to a |sin/c| kind of function. Data points in fc=0 or fc = tt diverges 
from the above sinusoidal function. This may be a finite size effect. In the thermodynamic limit 
of infinite chain length excitation spectrum of ferric wheel will be similar to Eq. ( pT| ) . 

We have also calculated the spin-spin correlation function {J2ij < ^i^j >) of ferric wheels and 
Fourier transformed it to find out the structure factor {S{q)). In Fig 5. (a - d) we have plotted the 
structure factor as a function of wave vector q for diff'erent symmetry subspaces. In all the cases 
S{q) shows a peak at g — tt. The ground state is a Neel ordered state. It also signifies that 
ground state will undergo a spontaneous distortion with wave vector tt. We find that the weight 
of S{q) is more in case of ground state than excited states. 
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B. Evolution of Magnetization in presence of an ac field 



We have studied the evolution of magnetization in presence of an axial ac magnetic field varying 
the amplitude of the field. We have kept the frequency of the field fixed at lo^IQ^^ . We calculate 
the magnetization at each time step. When we draw a smooth curve for the time evolution over 
long time periods, we find a sinusoidal motion, 

M{t) - cos{m) (12) 

which can be seen in Fig. 6. Unexpectedly the frequency Vl of this sinusoidal motion does not 
correspond to an eigenfrequency of the system or to the period of the external field. When we 
change the amplitude Hq of the field, the period of the magnetization changes, which is shown in 
Fig 6. We find a regular dependence of as a function of Hq, as in Fig 7. Frequency of oscillation Vl 
becomes very small for some values of Hq, again it increases but to a lower value than the previous 
^max- ^max valucs cau bc fitted to a exponentially decaying curve. This nontrivial resonance can 
be understood from the viewpoint of Floquet's theorenB. Nonadiabatic transitions are possible 
whenever system becomes metastable. This resonance can be observed if two states with different 
magnetizations are nearly degenerate leading to large tunneling cross sections between the states 
resulting in oscillation of the magnetization. This magnetization oscillation can be related to that 
of macroscopic quantum coherence, which is predicted by Zvezdin et. a/.lli. 



C. Torque Magnetometry 



Cornia et. al. have used a novel cantilever torque magnetometry technique to study the spin- 
state crossover in ferric wheels. The torque T experienced by a magnetically anisotropic substance 
in a uniform magnetic field B is given by 

T = M X B (13) 

where M is the magnetization of the sample. T vanishes when the magnetic field is applied along 
one of the principle directions x, y, z of the susceptibility tensor, since in this case M and B 
are collinear. The y component of the torque operator Ty can be easily obtained for an applied 
magnetic field of the form 

B = B(cos 6lz + sin 6lf), (14) 

in Eq. 5, and is given by, 

< Ty >= -gnBB{< Sx > COS0- < Sz> sin 6*) (15) 
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where < Sa > = '}Zi=i < ^i,a > IS thc ground state expectation value of the component a of the 
total spin operator. We have computed < T^^ > for Fei2 on the basis of the eigenvectors of the 
Hamiltonian in Eq. 5. When is j then torque is maximal and < Ty > is just the difference 
between the x and z spin components. In Fig. 8 we have shown the variation of < T,^ > with 
magnetic field. We can clearly observe the step behavior of the torque component which is a 
manifestation of the level crossing of singlet, triplet and quintet states. 



IV. SUMMARY AND OUTLOOK 



We have implemented a general and efficient procedure that allows us to block factorize the spin 
Hamiltonian matrix based on its invariance under cyclic symmetry and parity operation. This 
method can be used in general for systems of other symmetries also. We have obtained the low- 
lying eigenvalue spectrum of ferric wheels up to Fei2 using the C„ rotational symmetry of the 
molecules. We have also analyzed the dispersion spectrum and structure factor. To reproduce the 
low temperature properties of ferric wheels we need to know the low-lying eigenvalue spectrum of 
these systems. We have also studied the dynamics of ferric wheel evolving a proper initial state. 
We observe nontrivial oscillation of magnetization in presence of an alternating magnetic field. We 
have also obtained the torque of the ferric wheels in the presence of a nonaxial magnetic field and 
find that the torque also exhibits step like behavior with field. Evidently a study including the 
effect of nonzero temperature on this oscillation is a challenging problem of future research. 
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Figure Captions 



1. Representative Mg = state in (a) 6 spin-i cluster (b) Fei2 wheel with all the sites having spin 
S = |. Numbers in paranthesis correspond to the Mg value at the site. Bit representation as well 
as the integer value are given just below the diagrams. 

2. Plot of ground state energy (in unit of J) vs. inverse system size for ferric wheels up to Fei2. 

3. Plot of lowest energy (in unit of J) in every total spin sector vs. the total spin (Sjot) in case of 
Fei2. It clearly obeys Taft's conjecture. 

4. Plot of energy (in unit of J) vs. momentum vector k for Fcg, Fcg, Fcio and Fci2. 

5. (a) Plot of static structure factor {Sq (in arbitary units)) for Fee in all the subspaces. Here 
ground state lies in B+ subspace. In (b), (c) and (d) we have plotted the static structure factor 
in A and B subspaces for Feg, Feio and Fei2 respectively. In the doubly degenerate subspaces the 
weight of S(q) is much smaller than that in the ground state. 

6. Plot of evolution of magnetization in presence of an alternating axial magnetic field of three 
different amplitudes (Ho(in unit of 3/fi)). 

7. Field amplitude (Hq) dependence of frequency O (defined in the text). Dotted line shows the 
exponential nature of O^ax- 

8. Plot of y component of torque (in unit of ^_'^^^ )with change in magnetic field applied at ^ = 7r/4. 
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Energy gap vs. inverse system size 
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Fig. 2 
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Fig. 3 
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Fig. 5 (b) 
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Fe^o Ring 
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Fig. 5 (c) 
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Fig. 6 
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